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ABSTRACT 

We report on the discovery of thermal and non-thermal X-rays from the shells of the super 
bubble (SB) 30 Dor C in the Large Magellanic Cloud (LMC). The X-ray morphology is a nearly 
circular shell with a radius of ^40 pc, which is bright on the northern and western sides. The 
spectra of the shells are different from region to region. The southern shell shows clear emission 
lines, and is well fitted with a model of a thin-thermal plasma (kT — 0.21 keV) in non-equilibrium 
ionization (NEI) plus a power-law component. This thermal plasma is located inside of the Ha 
emission, which is the outer edge of the shell of the SB. The northern and western sides of the SB 
are dim in Ha emission, but are bright in non-thermal (power-law) X-rays with a photon index 
of 2.1-2.9. The non-thermal X-ray shell traces the outer boundary of the radio shell. These 
features of thin-thermal and non-thermal X-rays are similar to those of SN 1006, a prototype of 
synchrotron X-ray shell, but the non-thermal component of 30 Dor C is about ten-times brighter 
than that of SN 1006. 30 Dor C is the first candidate of an extragalactic SB, in which energetic 
electrons are accelerating in the shell. The age is much older than that of SN 1006, and hence 
the particle acceleration time in this SB may be longer than those in normal shell-like SNRs. 
We found point-like sources associated with some of tight star clusters. The X-ray luminosity 
and spectrum are consistent with those of young clusters of massive stars. Point-like sources 
with non-thermal spectra are also found in the SB. These may be background objects (AGNs) 
or stellar remnants (neutron stars or black holes). 

Subject headings: acceleration of particles — ISM: individual (30 Dor C) — X-rays: ISM 



1. Introduction 

Since the discovery of cosmic rays (Hess 1912), 
the sites and mechanisms of cosmic ray acceler- 
ation up to the knee energy, which is a spectral 
kink of cosmic rays at ~ 10 15 eV, have been key 
problems. The most plausible scenario is diffu- 
sive shock acceleration (DSA) in the shock fronts 
of supernova remnants (SNRs) (Bell 1978; Bland- 
ford & Ostriker 1978); in this scenario, particles 
receive energy whenever they cross the shock from 
downstream to upstream. 

Koyama et al. (1995) discovered synchrotron X- 
rays from the shock front of SN 1006, and con- 
firmed that the SNR accelerates cosmic rays up 



to ~TeV, together with the discovery of inverse 
Compton TeV 7-rays (Tanimori et al. 1998), which 
helps us to estimate the magnetic field and the 
maximum energy of electrons. Non-thermal X- 
rays are the best tracer of cosmic ray accelera- 
tion sites. In fact, several Galactic SNRs have 
already been confirmed as cosmic ray accelerators 
using hard X-rays (SN 1006; Koyama et al. 1995, 
Long et al. 2003, Bamba et al. 2003b, G347.3-0.5; 
Koyama et al. 1997, Slane et al. 1999, RCW 86; 
Bamba et al. 2000, Borkowski et al. 2001b, Rho 
et al. 2002, G266.2-1.2; Slane et al. 2001, Cas A; 
Vink & Laming 2003, and Tycho; Hwang et al. 
2002). The ASCA Galactic plane survey has also 
found some candidates (Bamba et al. 2001; Ueno 
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ct al. 2003; Bamba et al. 2003a). Although sev- 
eral Galactic SNRs have been confirmed to be sites 
of cosmic ray acceleration, some problems remain. 
For example, no SNR has been found in which cos- 
mic rays are accelerated up to the knee energy. In 
fact, Reynolds & Keohane (1999) have suggested 
that the maximum energy of cosmic rays acceler- 
ated in SNRs is at most ~ 10 14 eV, an order of 
magnitude below the knee energy. 

Super bubbles (SBs) are large hot cavities cre- 
ated by the combined actions of fast stellar winds 
and successive supernova (SN) explosions of mas- 
sive stars in OB star associations. Thus, SBs 
are strong candidates for being sites where cos- 
mic rays are accelerated to higher energies (c.f., 
Bykov & Fleishman 1992, Klepach, Ptuskin, & 
Zirakashvili 2000). However, it is difficult to 
study SBs in detail. Their distances, and thus 
their sizes and energies, remain uncertain in our 
Galaxy. Furthermore, interstellar absorption often 
prevents observations at optical, UV, and X-ray 
wavelengths. The Large Magellanic Cloud (LMC) 
provides an ideal location to examine the acceler- 
ation of cosmic rays in SBs because of its prox- 
imity (50 kpc; Feast 1999), low inclination angle 
(30°-40°; Westerlund 1997), and low foreground 
absorption (Ay < 0.3 mag; Bessel 1991). 

The LMC provides a sample of numerous SBs 
at a common distance that is resolvable by mod- 
ern X-ray detectors. For this study, we singled out 
30 Dor C as being the best target for our search. 
A radio source was discovered by Le Marne (1968) 
south-west of 30 Dor, and was named 30 Dor C, 
which is now categorized as a SB. Mills ct al. 
(1984) found a shell-like structure in the 843 MHz 
band observation with a radius of about 3' (—40 pc 
on 50 kpc distance). Along the radio shell, com- 
plex Ha, HP, and [SII] emissions were found, 
which are bright in the south-east, but dim in the 
other side (Mathewson et al. 1985). In the X-ray 
band, Einstein detected the eastern shell of the 
30 Dor region for the first time (Long, Helfand, & 
Grabelsky 1981; Chu & Mac Low 1990). Dunne, 
Points, & Chu (2001) (hereafter, DPC) reported 
that the ROSAT spectrum required a thermal 
model with a rather high temperature (~1 keV). 
Itoh et al. (2001) found non-thermal X-rays in the 
ASCA data, but could not spatially resolve the 
thermal and non-thermal components in the SB. 
Recently, Dennerl et al. (2000) found a complete 



ring in soft and hard X-rays with a diameter of 
~6'with XMM Newton. The hard X-rays from the 
shell resemble the synchrotron X-rays in the shock 
front of SNRs, which is a site of cosmic ray acceler- 
ation. Therefore, 30 Dor C would be the first and 
good candidate of cosmic ray accelerating SBs. 

In this paper, we report on the first results of 
spatially resolved hard X-ray spectral analyses of 
30 Dor C using Chandra and XMM-Newton data. 
In §2, we describe the data and their reduction. 
We analyze the data in §3, followed by a discussion 
in §4, and a summary in §5. The distance to the 
LMC is assumed to be 50 kpc (Feast 1999). 

2. Observations and Data Reduction 

Chandra observed SN 1987A several times with 
the Advanced CCD Imaging Spectrometer (ACIS) 
S array; the satellite has an excellent spatial res- 
olution (~0.5" on the aim-point), and one CCD 
has a field of view (FOV) of ~ 8.5' x 8.5', an en- 
ergy range of 0.3-10.0 keV, and an energy res- 
olution of -130 cV at 5.9 keV. The details for 
the satellite and the instruments are described by 
Garmire et al. (2000) and Weisskopf et al. (2002), 
respectively. Some of the SN 1987A observations 
covered the 30 Dor C region by chance, which is 
located at about 5' northeast from SN 1987A. In 
order to study the diffuse structure of 30 Dor C, 
we selected two observations in which the grat- 
ing instrument is removed from the X-ray mirror 
axis (Observation ID = 1044 and 1967; hereafter 
Obs.l and 2). The observed dates and targeted 
positions are given in Table 1. Data acquisition 
from the ACIS was made in the Timed-Exposure 
Faint mode with a readout time of 3.24 s in both 
observations. The data reductions and analyses 
were made using the Chandra Interactive Analysis 
of Observations (CIAO) software version 2.3. For 
a Level- 1 processed events provided by a pipeline 
processing at the Chandra X-ray Center, we made 
a Charge Transfer Inefficiency (CTI) correction, 
and selected ASCA grades 0, 2, 3, 4, and 6, as the 
X-ray events. High-energy events due to charged 
particles and hot and flickering pixels were re- 
moved. There are some "streaks" in the CCD 
chip S4, which are caused by a random deposition 
of a significant amount of charges in pixels along 
the row during the serial read-out process. These 
"streak" events were removed by using the pro- 
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gram destreak 1 in CIAO. The total available times 
for each observation, after the screening, are listed 
in Table 1. 

XMM-Newton (Jansen et al. 2001) also ob- 
served SN 1987A and the 30 Dor region several 
times; the satellite has a spatial resolution of 
<~14" and a relatively wide FOV with the radius 
of 15'. We selected two observations (Observa- 
tion ID = 0104660301 and 0113020201; hereafter 
Obs.3 and 4), which cover 30 Dor C and are rel- 
atively free from high background flares due to 
low-energy protons. The observed dates and tar- 
geted positions are shown in Table 1. In both 
observations, only the metal oxide semiconduc- 
tors (MOS) CCDs, which have an energy range 
of 0.1-10.0 keV and a similar energy resolution to 
Chandra, (Turner et al. 2001) were operated in the 
full-frame mode with the medium filter (Stcphan 
ct al. 1996; Villa et al. 1998) for blocking ultra- 
violet photons. The data reductions and analy- 
ses were made using the XMM-Newton Standard 
Analysis System (SAS; Watson ct al 2001) ver- 
sion 5.4.1; we performed the basic pipeline pro- 
cess following the SAS guide. The background 
level was largely changed, particularly in Obs. 4, 
hence we removed the data with a high back- 
ground level (>0.6 cnts s _1 in the 10.0-15.0 keV 
band). The exposure times in each observation 
after the screenings are listed in Table 1. 

3. Analyses and Results 

3.1. Diffuse Emission and Point-like Sources 

Figure 1 shows the soft (0.7-2.0 keV) and 
hard (2.0-7.0 keV) band Chandra images around 
30 Dor C, in which the two observations (Obs.l 
and 2) arc combined with a correction of the ex- 
posure times. A clear shell-like structure with a 
radius of ~ 170" (R <~ 40 pc-radius at the 50 kpc 
distance) is seen in both bands. In detail, however, 
the morphologies are different from each other; 
the entire shell is seen in the soft band, whereas 
the hard X-rays are visible only at the western 
part. Catalogued SNRs, the Honeycomb nebula 
(SNR 0536-69.3) and SN 1987A are also seen 
mainly in the soft band (see Figure 1). The X-ray 
features of these objects have been reported with 
XMM-Newton (Dennerl et al. 2000) and Chandra 

See http://asc.harvard.edU/ciao2.3/ahelp/destreak.html. 



(Burrows et al. 2000; Park et al. 2002; Michael et 
al. 2002) monitoring observations. 

In addition to the diffuse structure, some point- 
like sources are found inside 30 Dor C. The expo- 
sure time of the XMM-Newton observations are 
short, most of the data suffer from a high back- 
ground, and the spatial resolution is not suffi- 
cient. Hence, we concentrated on the Chandra 
data for the point-source search and analysis. At 
first, point sources were searched for with the 
wavdetect softwarein the 0.5-8.0 keV band images, 
then manually inspected for any spurious point- 
like structure due mainly to a part of the diffuse 
emission. We thus found six point sources with 
a significance level of >7.0a, as shown in Fig- 
ure 1 and Table 2. For these six point sources, we 
searched for optical, infrared, and radio counter- 
parts, and found that three (No.l, 3, and 4) coin- 
cide at the positions of the brightest star clusters: 
a, /?, and 7 (Lortet & Testor 1984). We therefore 
checked further for any X-ray emission from the 
other clusters (5, e, and £), which are also mem- 
bers of the OB association LH 90 (Luckc & Hodge 
1970) encompassed by the 30 Dor C shell. How- 
ever, we found no excess X-rays from these clusters 
at the 3cr limit. We also searched for X-ray coun- 
terparts from ROSAT PSPC and HPJ catalogues 
(Haberl & Pietsch 1999; Sasaki, Haberl, & Pietsch 
2000), but found no candidate. 

3.2. Spectra 

In the Chandra observations, we can see that 
the diffuse structure consists of several shell frag- 
ments (see Figure 1), and the whole structure is 
widely spread over the two observed regions with 
different configurations of the CCD types (back- 
illuminated and front-illuminated). Therefore, a 
spectral analysis on all of the diffuse structure is 
technically and scientifically complicated. For this 
reason, we divided the diffuse structure into four 
regions (hereafter shells A-D), as shown in Fig- 
ure 1, and performed the spectral analysis sepa- 
rately for each shell. We excluded all of the de- 
tected point sources (see §3.1) from each shell. 
The background regions were also selected sepa- 
rately for each shell from the source-free regions 
in the same CCD and observations as those of the 
shells. 

The XMM-Newton fields covered the whole dif- 
fuse structure. However, for consistency with the 
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Chandra analysis, we divided the diffuse struc- 
ture into shells A-D, the same as in the Chan- 
dra case. The background regions were selected 
from the source-free regions in the same observa- 
tion. The spectra of each shell taken from all the 
available observations with Chandra and XMM- 
Newton were analyzed simultaneously. However, 
for brevity, we show only the Chandra spectra 
(background-subtracted) and the fitting results in 
Figure 2. 

For spectral analyses, we used XSPEC vll.00 
(Arnaud 1996) The spectrum of shell A shows 
many line-like structures with the center energies 
at 0.58, 0.68, 0.92, and 1.35 keV, which correspond 
to the emission lines of He-like O Ka and K/3, He- 
like Ne, and He-like Mg, respectively. We there- 
fore fit the spectrum with a thin-thermal plasma 
model in non-equilibrium ionization (NEI) (NEI; 
Borkowski et al. 2001a) with the mean LMC abun- 
dances (Russel & Dopita 1992; Hughes, Hayashi, 
& Koyama 1998). The interstellar absorption in 
our Galaxy and LMC were treated separately. 
The Galactic absorption column was estimated us- 
ing the HI data by Dickey & Lockman (1990) as 
-Nh.hi = 6.35 x 10 20 cm~ 2 . Arabadjis & Bregman 
(1999) reported that the value of Na, measured in 
the X-ray band, is twice that of A^h.hi in the case 
of \b\ > 25° and 7V h ,hi > 5 x 10 20 cm" 2 . There- 
fore, we fixed the galactic absorption column to 
be A^h = 1.27 x 10 21 cm~ 2 ; we used the cross sec- 
tions by Morrison & McCammon (1983) and the 
solar abundances (Anders & Grevesse 1989). The 
absorption column in the LMC was, on the other 
hand, treated as a free parameter with the mean 
LMC abundance (Russel & Dopita 1992; Hughes, 
Hayashi, & Koyama 1998). 

This thin-thermal plasma model was rejected 
with a x 2 /degree of freedom (d.o.f.) = 331.8/211, 
even if we allow the abundances to be free, leav- 
ing a systematic data residual at the high-energy 
band. We hence added a power-law component on 
the thin-thermal model. Since the two-component 
model still leaves large residuals at about 0.8 keV 
and 1.3 keV, we allowed the abundance of Fe and 
Mg in the thermal plasma ( "NEI" component) to 
be free. The fitting was then greatly improved 
with a x 2 /d.o.f = 264.6/209. Although this two- 
component model is still rejected in a statistical 
point of view, further fine tuning on the model 
is beyond the scope of this paper. Figure 2 and 



Table 3 show the best-fit models and parameters, 
respectively. 

Unlike shell A, the X-ray spectra of shells B- 
D are hard and featureless, suggesting non- 
thermal origin. In fact, a thin thermal model 
fitting requires an unrealistically high tempera- 
ture (> 2 keV) and low abundances (z < 0.3). 
We therefore fitted the spectra with a power-law 
model with absorption, which was calculated in 
the same way as that for shell A, and found ac- 
ceptable fits for all of the spectra. The best-fit 
models and parameters are shown in Figure 2 and 
Table 4, respectively. 

It is conceivable that the spectra of shells 
B-D may include a small fraction of the thin- 
thermal component found in shell A. We therefore 
added the same thin-thermal spectrum as that 
for shell A, and fitted with this composite model 
(thin-thermal plus power-law). The free param- 
eters are normalizations of the two components: 
power-law index and Nn value. However, no sig- 
nificant fraction of the thin thermal component is 
found from shells B-D. 

For all of the point sources, the X-ray photons 
are collected from an ellipse with the radii of the 
point spread function (PSF), as listed in Table 2. 
We note that all of the sources are located far from 
the on-axis position of the X-ray mirror, and hence 
the PSFs are larger than the best value of Chandra 
PSF (~ 0.5"on the aim point). The background 
regions were selected from source-free regions in 
the same way as the diffuse emissions. We first fit- 
ted the spectra with a thin-thermal plasma model 
in collisional equilibrium (MEKAL; Mewe, Gronen- 
schild, & van den Oord 1985; Kaastra 1992) with 
an absorption calculated in the same way as diffuse 
emission. The abundances are fixed to be 0.3 so- 
lar, the average value of interstellar medium in the 
LMC. The fittings are acceptable for two sources 
(No. 1 and 4) with reasonable temperature (2.1 
and 1.0 keV), but for the spectra of the other 4 
sources, the models are either rejected or require 
an unreasonably high temperature. We therefore 
fitted the spectra of these sources with a power-law 
model. The best-fit parameters and the reduced 
X 2 s are listed in Table 2. 
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4. Discussion 

4.1. The absorption 

The absorption columns of the north-eastern 
shells (shells A and B) are similar to those of 
most sources in the LMC (~ 10 21 cm -2 ), whereas 
those of the other shells (C and D) are signifi- 
cantly larger (~ 10 22 cm -2 ) than the typical LMC 
absorption. A similar trend was found for the 
point sources; those in the western half, No.l and 
2, have a larger (~ 10 22 cm -2 ) than those 
in the eastern sources, No. 3-6 (< 10 21 cm -2 ). 
Since DPC (2001) has already reported this ten- 
dency, we have thus confirmed the results with the 
better spatial and spectral capability of Chandra. 
This systematic increase of absorption toward the 
western region of 30 Dor C may be due to extra 
absorption of a molecular cloud located in front 
of the western half of 30 Dor C. To verify our 
conjecture, we searched for the molecular cloud 
in the CO map (Figure 2(a) in Yamaguchi et al. 
2001) and found the candidate with the intensity 
of /(CO) ~ 3.6 K km s -1 . With a conversion fac- 
tor of AT(H 2 )/J(CO) ~ 9xl0 20 cm -2 (K kms -1 ) -1 
(Fukui et al. 1999), the estimated absorption col- 
umn due to the molecular cloud is N^ IC ~ 6.5 x 
10 21 cm -2 , which is consistent with our result. 

4.2. The Thermal Emission 

30 Dor C is a SB made by a strong stellar 
wind and/or successive supernova explosions of 
massive stars located in the OB star association 
LH 90. The age of this star association, or that 
of 30 Dor C, is on the order of a few to 10 Myr 
(Luckc & Hodge 1970). The thermal emissions are 
enhanced in the south-eastern side of 30 Dor C 
(around shell A). The position of this compo- 
nent coincides with that of the Ha emission (DPC 
2001). The plasma temperature of 30 Dor C 
is rather high compared to those of other LMC 
SBs (DPC 2001), although it becomes significantly 
lower than previous results (DPC 2001). Perhaps 
due to the poor spectral resolution of ROSAT, 
previous observations could not resolve the power- 
law component (hard spectrum) from the thermal 
emission. 

The X-ray luminosity of 30 Dor C is signifi- 
cantly lower than that of the other SBs (DPC 
2001). With the assumption that the plasma 



in shell A distributes uniformly in the ellipsoid 
with radii of 58" x 34" x 34" (total volume V = 
1.2 x 10 59 cm 3 ), the mean density (n c ), thermal 
energy (E) and the age of the plasma (t p ) were 
calculated as follows using the emission measure 
E.M. = n^V and ionization time scale (see Ta- 
ble 3): 

n c = 4.7 (4.2-4.8) x 10 -1 [cm -3 ], (1) 
E ~ 3n c kTV = 5.7 (4.6 - 6.4) x 10 49 [erga]2) 
i p = 6.7 (> 0.9) x 10 5 [yrs]. (3) 

Together with the Ha emission around shell A 
(DPC 2001), we infer that the shell of 30 Dor C 
collides with dense matter and temporally emits 
thermal X-rays with a relatively high temperature. 
The overabundance of a light element (Mg) rela- 
tive to heavier element (Fe) may indicate that a 
type II SN occurred (Tsujimoto et al. 1995). Thus, 
the progenitor is a massive star which is a member 
of cluster LH 90 (Luckc & Hodge 1970) near the 
center of 30 Dor C (see section §4.4). 

4.3. The Non-Thermal Emission 

In the 843 MHz band, Mills et al. (1984) found 
a clear radio shell, which is brightest on the south- 
western side (around shells C and D) and dim on 
the eastern side. The non-thermal X-ray emissions 
are enhanced at the radio bright shell. This fact 
implies that the non-thermal X-rays are emitted 
by the same mechanism as the radio band. There- 
fore, the non-thermal X-rays are likely to be syn- 
chrotron radiation from the accelerated electrons, 
like SN 1006 (Koyama et al. 1995) and other SNRs, 
which accelerate particles up to TeV. The X-ray 
photon index of 2.1-2.9 (see Table 4) is, in fact, 
typical to synchrotron emissions. 

To verify the synchrotron origin, we fitted the 
X-ray spectra with a SRCUT model, which repre- 
sents the synchrotron emission from electrons with 
an energy distribution of a power-law plus expo- 
nential cut-off (Reynolds 1998; Reynolds & Keo- 
hane 1999). Since the radio index data is not accu- 
rate due to the large background and contamina- 
tion of thermal emission (Mathewson et al. 1985), 
we fixed the spectral index (a) at 1 GHz to be 0.5, 
which is expected from the first-order Fermi accel- 
eration and similar to that of SN 1006 (a = 0.57; 
Allen, Petre, Gotthelf 2001). 

The fittings were statistically acceptable, and 
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the best-fit parameters are listed in Table 4. The 
best-fit cut-off frequency is also similar to that 
of SN 1006, although the age of 30 Dor C may 
be on the same order as that of LH 90 (10 Myr; 
Lucke & Hodge 1970), which is far larger than 
SN 1006. This implies that, depending on the 
environment, the acceleration of high-energy elec- 
trons up to the knee energy can continue for a 
far longer time than the previous consensus for 
the SNR case (~ 10 3 yrs; Reynolds and Keohane 
1999). The electron acceleration time may have 
been extended, because successive supernova ex- 
plosions in 30 Dor C, possibly over the course of a 
few Myr, may more-or-less continuously produce 
high-energy electrons. The total luminosity of the 
non-thermal component (<~ 5.3 x 10 35 ergs s _1 ) 
is about 10-times larger than that of SN 1006 
(Koyama et al. 1995). This large non-thermal flux 
in 30 Dor C would be due to the large energy sup- 
ply by multiple successive supernova explosions. 

The expected flux density at 1 GHz is 3.0 (2.4- 
4.6) xlO -2 Jy with a = 0.5, which is significantly 
smaller than the observed value (1.0 Jy; Mills et al. 
1984). This "inconsistency" is not relaxed, even if 
we assume a larger radio index of a = 0.6 (1.8 
(1.4-2.5) xlO^Jy, as shown in Table 4). We in- 
fer that the larger observed radio flux than that 
expected from X-rays would be due to either the 
contamination of thermal radio flux or a large un- 
certainty of the background level (see Mills et al. 
1984). 

The Ha emission (DPC 2001) is anti-correlated 
with the non-thermal components. Similar fea- 
tures have been observed in some of other SNRs 
with synchrotron X-rays: SN 1006 (Winkler & 
Long 1997) and RCW 86 (Smith 1997). Because 
the Ha region has higher density, it may have a 
higher magnetic field. Therefore, the maximum 
electron energy is limited by the quick synchrotron 
energy loss, leading to reduced non-thermal X- 
rays. 

4.4. Point Sources 

We have identified three point-like X-ray sources 
(No. 1, 3, and 4) to the tight clusters of massive 
stars a, (3, and 7, respectively. The spectrum 
of No.l (a) is fitted with a thin thermal plasma 
model of ^2.1 keV, which is consistent with stel- 
lar X-rays from massive stars. The cluster a is 
the brightest of the three, with X-ray luminos- 



ity of ~ 10 34 ergs s 1 . Optical spectroscopy of 
this cluster revealed that it includes one red gi- 
ant and one Wolf-Rayet (WR) star (MG 41 and 
Brey 58; Lortet & Testor 1984). Hence, the latter 
is a possible counterpart of the 2 keV source. The 
X-ray luminosity is near to the upper end of a 
massive star, or its binary (Maeda et al. 1999). 
Although the X-ray luminosity of No. 3, a counter- 
part of the star cluster (3, is also consistent with 
a massive young star, the best-fit photon index, 
^2 (see Table 2), is rather typical to a rotation 
powered neutron star, probably a stellar remnant 
of an SN explosion in the active star cluster. The 
spectrum of No. 3 (7) is soft and well-fitted with 
a, kT = 1.0 keV thermal plasma model, similar 
to a cluster a. The star cluster 7 also includes 
one OB star (Sk -69?212; Sanduleak 1970), and 
the X-ray luminosity is consistent with that of a 
young massive star. 

The other three sources (No, 2, 5, 6) have 
no counterpart in optical, in the infrared band 
(Breysacher 1981; Lortet & Testor 1984) nor in the 
SIMBAD data base. Their spectra are relatively 
hard, and are consistent with being background 
active galactic nuclei (AGNs), or stellar remnants 
(black hole or neutron star) by successive SN ex- 
plosions. 

In order to constrain the nature of the point- 
like X-ray sources with a power-law spectrum, 
we further examined the time variability with 
the Kolmogorov-Smirnov test (Press et al. 1992). 
However, no significant time variability was found 
even between the two observations. For high reso- 
lution timing, we examined the Chandra High Res- 
olution Camera (HRC) data (ObsID = 738), but 
no pulsation was found from any of these point 
sources. 

5. Summary 

(1) In the shell of the super bubble 30 Dor C, 
we resolved non-thermal and thermal X-rays using 
Chandra and XMM-Newton data. 

(2) The thermal emission concentrates on 
the south-eastern shell of 30 Dor C. The spec- 
trum is well represented by a thin thermal 
plasma model with kT =0.21 (0.19-0.23) keV and 
n c =0.47 (0.42-0.48) cm" 3 . 

(3) The non-thermal X-rays are located at 
northern and western parts of the SB. The power- 
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law model is well-fitted to the spectra with 
r =2.1-2.9, similar to that of the cosmic ray accel- 
erating SNR, SN 1006. This is the first discovery 
of non-thermal X-rays from the shells of SBs. The 
total luminosity is ten-times larger than that of 
SN 1006. 

(4) We found six point-like sources in 30 Dor C. 
Three sources are located in tight clusters of mas- 
sive stars, and therefore they may be active Wolf- 
Rayet stars or compact stars. The other three 
have power-law spectra with no optical counter- 
part, implying that they are background AGNs or 
compact remnants of SNc. 
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Fig. 1.— Chandra ACIS images around 30 Dor C in the 0.7-2.0 keV band (a) and in the 2.0-7.0 keV band 
(b), with J2000 coordinates. The scales are logarithmic with the units of xl0~ 5 counts s _1 arcmin _2 cm -2 , 
as shown in the left bar for each image. The source regions used for the spectral analyses are shown by the 
solid lines in the left panel (Fig. 2, Table 3 and 4). The crosses represent the position of the detected point 
sources shown in Table 2. The position of SN 1987A (Burrows ct al. 2000; Park et al. 2002; Michael et al. 
2002) and the Honeycomb nebula (Dcnncrl ct al. 2000) are also shown. 



Table 1 

Observation log around 30 Dor C. 





Satellite 


Obs ID 


Position (J2000) 




Date 


Exposure 








(RA, DEC) 




(yyyy/mm/dd) 


(ks) 


Obs.l 


Chandra 


1044 


(05 h 32 m 22 i :0, -69 d 16" 


33?0) 


2001/04/25 


18 


Obs.2 


Chandra 


1967 


(05 h 32 m 22!0, -69 d 16" 


'33 B 0) 


2000/12/07 


99 


Obs.3 


XMM-Newton 


0104660301 


(05 h 35 m 28?0, -69 d 16" 


11?0) 


2000/11/25 


21 


Obs.4 


XMM-Newton 


0113020201 


(05 h 37 m 47!6, -69 d 10" 


20=0) 


2001/11/19 


16 
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Shell A - 







I 7 o 



channel energy (keV) 





PV. 


Shell B - 













5 1 


2 5 
channel energy (keV) 


10 



Shell C 





± 




channel energy (keV) 



channel energy (keV) 



Fig. 2. Upper panels: Background-subtracted spectra for the shells A-D (crosses). Although the spectral 
fittings were made with both the Chandra and XMM-Newton data simultaneously, only the Chandra data 
and results are shown for brevity. The best-fit models are shown with solid (power-law) and dotted (NEI) 
lines. Lower panels: Data residuals from the best-fit models. 



Table 2 

Point source data around 30 Dor C. a 



No. 


Name 
CXO JOS... 


PSF b 


SM C 


F/kT 
■■■ /[keV] 


TirLMC d 

[xlO 22 Hem -2 ] 


Flux' 
[ergs cm 


-in 


Luminosity 
Teres s _1 l 


X 2 /d.o.f. 


1 . . . 


3542.4-691152 


7.0 x 3.6 


MK 


2.1 (1.4-3.6) 


3.5 (2.5-5.1) 


1.8 X 10" 


-14 


1.2 x 10 34 


9.53/12 


2 ... 


3542.9-691206 


7.3 x 3.8 


PL 


2.1 (1.4-3.2) 


7.7 (4.4-13) 


2.0 X 10" 


-14 


1.5 xlO 34 


8.25/13 


3 ... 


3559.9-691150 


7.1 x 4.5 


PL 


1.7 (1.3-2.6) 


0.13 (< 1.3) 


6.5 x 10" 


-15 


2.2 xlO 33 


4.03/7 


4 . . . 


3606.6-691147 


8.3 x 6.6 


MK 


1.0 (0.19-1.3) 


(< 1.9) 


2.8 x 10" 


-15 


1.1 xlO 33 


8.64/9 


5 . . . 


3620.7-691303 


8.5 x 4.4 


PL 


1.9 (1.5-2.4) 


(< 0.22) 


8.9x10" 


15 


3.1 xlO 33 


19.8/20 


6 ... 


3633.3-691140 


10.7 x 8.5 


PL 


1.8 (1.5-2.3) 


0.5 (0.08-1.1) 


2.0 x 10" 


-14 


7.8 xlO 33 


21.4/18 



a Parcnthcscs indicate single-parameter 90% confidence regions. 
b Thc elliptic PSF size. 

c Spcctral model; PL: power-law, MK: thin thermal plasma in collisional equilibrium (Mcwc, Groncnschild, & van den Oord 1985; 
Kaastra 1992) with the abundance of 0.3 solar. 

d Absorption column in the LMC. The abundance is assumed to the average LMC values (Russel & Dopita 1992; Hughes, Hayashi, 
& Koyama 1998). 

°Absorbed flux in the 0.5-9.0 keV band. 

f Intrinsic luminosity (absorptions are removed) in the 0.5-9.0 keV band, at the LMC distance of 50 kpc (Feast 1999). 
NOTE. — No.l, 3 and 4 have the optical counterparts, the stellar clusters a, j3, and 7 respectively (Lortet & Testor 1984). 
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Table 3 

Best-fit parameters of the shell A. a 



Parameters NEI Power-law 



fcT/r [kcV/— ] 0.21 (0.19-0.23) 2.9 (2.7-3.1) 

n o i p b [10 12 cm~ 3 s] 9.9 (>1.4) 

_EM C [10 58 cm -3 ] 2.6(2.1-2.8) 

[Mg/H] d 3.3 (2.5-4.0) 

[Fe/H] d . 0.10 (0.06-0.14) 

Flux°[ergs cm^s -1 ] 6.8 X 10~ 14 8.4 X 10~ 14 

Af^ MCt [xl0 22 cm" 2 ] 0.24 (0.21-0.28) •■• 8 



a Parentheses indicate single-parameter 90% confidence regions. 

b Ionization time-scale, where n c and t p are the electron density 
and age of the plasma. 

c Emission measure EM — f n c n p eZV ~ n^V ~ n c^i where n p 
and V are the proton density and the plasma volume, respectively. 
The distance to LMC is assumed to be 50 kpc. 

d Abundancc ratio relative to the solar value (Anders & Grevesse 
1989). 

"Absorbed flux in the 0.5-9.0 keV band. 

f Absorption column in the LMC. The abundances are assumed 
to be the average LMC values (Russcl & Dopita 1992; Hughes, 
Hayashi, & Koyama 1998). 

8 Fixed to the same value as that for the NEI component. 



Table 4 

Best-fit parameters of the shells B-D. 



b 



Power-law 

Photon Index 2.7 (2.5-2.9) 2.3 (2.1-2.4) 2.5 (2.3-2.7) 

Flux^ergs s^'em" 2 ] 4.7 Xl0~ 13 4.5 Xl0~ 13 1.5 Xl0~ 13 

JV™ C [10 22 cm~ 2 ] c . 0.10 (0.06-0.19) 1.1 (0.94-1.4) 1.0 (0.85-1.3) 

X 2 /d.o.f 193.0/138 68.3/84 121.5/110 

SRCUT (a = 0.5) 

^roiioff [10 16 Hz] 6.2 (4.1-7.9) 24 (15-47) 12 (6.4-23) 

Flux density at 1 GHz [10~ 2 Jy] 1.9 (1.5-3.0) 0.71 (0.67-0.75) 0.42 (0.25-0.83) 

N^ MC [10 22 cm -2 ] . (<0.04) 1.0 (0.90-1.1) 0.89 (0.75-1.0) 

X 2 /d.o.f 195.5/138 67.0/84 120.2/110 

SRCUT (a = 0.6) 

^roiioff [10 16 Hz] 8.0 (5.3-11) 36 (20-75) 16 (8.3-35) 

Flux density at 1 GHz [10~ 2 Jy] 11 (8.2-16) 4.1 (3.9-4.3) 2.5 (1.6-4.6) 

JVh MO [10 22 cm _2 ] c . (<0.01) 1.0(0.91-1.2) 1.1(1.0-1.3) 

X 2 /d.o.f 195.2/138 67.1/84 121.8/110 



a Parentheses indicate single-parameter 90% confidence regions. 
b Flux in the 0.5-9.0 kcV band. 

c Absorption column in the LMC. The abundances arc assumed to be the average LMC values (Russcl 
& Dopita 1992; Hughes, Hayashi, & Koyama 1998). 
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